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In Classical Mechanics courses, the so-called Kepler problem pL s ch. 3, sec. 3.7] is 
virtually the only serious completely integrable application of central forces discussed 
besides the harmonic oscillator. As Sivarcliere already pointed out in this journal [2], the 
motion of a charged particle in the held of a magnetic monopole j3], which is another 
example of completely integrable problem, is, unfortunately, not discussed. 

Here, we extend and deepen Sivardiere’s study to a more general case, namely the 
interaction between Schwinger’s dyons (dual charged particles) [3]. As it has formal 
similarity with the Kepler problem and the Rutherford scattering, we believe that it 
may be presented right after these ones in Classical Mechanics courses. Furthermore, 
it exhibits unusual features such as non-planar stable orbits and rainbow and glory 
scattering, results that may arouse students’ and teachers’ interest. 

2. A brief historical review of magnetic monopoles and dyons 

The similarity between the electric and magnetic fields is visible in Maxwell’s equations. 
When deducing them in 1873, Maxwell himself pointed out that it would be necessary 
to assert that there are no net magnetic charged bodies and no ‘magnetic currents’ [5] 
art. 380, p. 6]. 

However, this symmetry would be restored if we were to assume the existence of a 
magnetic held density p m and a ‘magnetic current’ j m , obtaining (in SI units and with 
magnetic charges measured in ampere ■ meters) 


V • E 

Pe/f 0; 

V x E 

= -Pojm - dB/dt, 

V B 

hoPm; 

V x B 

= hoje + /i 0 eo dE/dt. 


In 1896, H. Poincare, applied the concept of magnetic matter to explain Birkeland’s 
magnetic deflection of cathode rays experiment. In it, he considered that the beam 
passes so close to one pole of the magnet that the other can be neglected [6]. This can 
be interpreted as a magnetic monopole approximation. 

Nevertheless, the idea of a magnetic monopole as a particle having a single magnetic 
pole was introduced only in 1931 by Dirac, in his famous work [3]. 

As a matter of fact, Dirac was not looking for something like the monopole, but 
investigating the why of the quantization of electric charge. In other words, why the 
electric charge always appears in Nature as a multiple of the electron charge e and why 
this charge has a value such that (in the same units as in ([Tj) ) 

2 e 0 hc/e 2 = 137. (2) 

In that work, however, instead of the relation ([2]), he obtained 

eg/e 0 hc 2 = Ji (n — ±1, ±2, ±3,...), (3) 
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which is known as the Dirac quantization condition. It does not set a value for e, but 
only for its product by the magnetic charge g of a hypothetical particle in its vicinity. 
On the other hand, Caruso |7j arrived at the same result (|3]) through a semi-classical 
derivation that gives it a new interpretation. 

Despite that, in a way, (J3|) offers a solution to Dirac’s initial problem: if there 
were a single magnetic monopole in the entire Universe, then all electric charges would 
be quantized according to ([3]). Therefore, in view of the observed quantization of 
electric charge and in the absence of another explanation for this fact at the time, 
the condition ([3]) was considered a serious argument for the existence of at least one 
magnetic monopole in the Universe. 

Generalizing Dirac quantization condition, Schwinger [S3], in 1969, introduced the 
dyons, which were suggested at the time as candidates for the quark model. Pinfold et 
al. [8J discuss the tremendous implications that the discovery of magnetic monopoles or 
dyons would have for our fundamental understanding of Nature at the deepest levels 
and describe the search for these particles, from Dirac’s proposal in 1931 to the current 
MoEDAL experiment at CERN’s LHC. 

After this short historical account, we proceed to study the interaction of two dyons 
from a classical mechanics point of view. 


3. The classical mechanics of two dyons 


To formulate the Lagrangian for the electromagnetic interaction of two dyons, we need to 
circumvent the problem of what became known as the Dirac string [3]. It is a singularity 
line in the space starting from the monopole, over which the magnetic vector potential A 
does not satisfy the condition that the wave function is a univalent function. 

Later on, Schwinger [9], Yang [10] , and Wu and Yang mi considered that the Dirac 
string has no physical meaning or real existence, being only an effect of the coordinate 
system, analogous to the problem of terrestrial geographic poles when trying to map 
the Earth’s surface with a single chart. These authors stated that this problem can be 
solved by dividing the space around the monopole in two regions a and b and defining 
two potential vectors A a and Af, which describe the monopole field in each of these 
regions and have singularities in the other b and a ones, respectively. 

Bollini and Giambiagi D2j. in their turn, proposed a multivalued distribution 
(generalized function) potential instead of a singular function. This approach, however, 
demands the use of the mathematical theory of distributions to the evaluation of its 
curls and divergences. 

On the other hand, Sokolov |T3] showed that the singularity of the magnetic 
monopole potential is of a purely kinematic origin, caused by the uncertainty of the 
azimuthal (ft angle along the z axis. As a consequence, for a charge in the field of a 
’Coulombian’ magnetic monopole 

n h-0 9 

B = - - -r 

47r r 6 


( 4 ) 
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we can obtain the equation of motion, by using the Lorentz force 
F = ev x B 


( 5 ) 


as 


and the expression (|4]) for the magnetic held, 
dv p 0 r 

» dt = ^ 69V X ^ 

without the appearance of strings or fictitious fields as long as vector potentials and 
curls are written and evaluated in spherical coordinates in 


( 6 ) 


dv .. 

h= ev x (V x A) spherical . 


(7) 


Sokolov’s procedure can be as well generalized to the dyon-dyon case ta p- is] 
for any vector potential A whose curl evaluated in spherical coordinates furnishes the 
correct magnetic potential Q such that ((Tj) leads to As a result, it is possible to 
obtain a classical string-free Lagrangian for the dyon-dyon case. 

That being said, we will here build the dyon-dyon Lagrangian through a different 
procedure, by means of the fields instead of a potential. This procedure will lead, 
however, to a Lagrangian that matches Sokolov’s one. 

We start by noticing that Maxwell’s equations ([Tj) are invariant under the duality 
transformation E —» cB, cB —> —E, p e —* p m /c, and j e —* j m /c. That allows us to 
generalize Lorentz force to the dyon-dyon interaction and write the equation of motion 
as 


dv . 
»Tt =e ' {E 


v x B) + gi ( B - -v x E 


= ei 


. ho 


r 

+ 


ho 

61 4^ 2 


1 1 


,9r 


47re ( 


-e 2 


r 

v x — 


47Te 0 

1 


eie 2 


r9i92 


r , ho f 
+ ( ' 192 


9ie 2 ) v x — 


"9 


47 T 60 r 3 47T 


ho r 

KV X — 


( 8 ) 


where 


(9) 


q = eie 2 + gig 2 /c 2 

K ~ e l92 ~ 9l e 2, 

being ei, e 2 , gi, and g 2 the electric and magnetic charges of the two dyons, corresponding 
the index 2 to the dyon that remains at the origin of the relative coordinate system, c 
the speed of light in a vacuum, and //, naturally, the system reduced mass, given by 

fflim 2 

h = -:-• 

mi + m 2 

To evaluate the cross product in 
coordinates as (see Figure [I]) 

f* — 

v — ff + r96 + r sin Q<p (j). 


( 10 ) 

we need to express r and v in spherical 


( 11 ) 
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Now, we can write the force on the right-hand side of (|8|) as 


F = 


i q 

47 T 60 r 2 

i q 

47ren r 2 


ho 

r - — 
47r 

r + &, 


K . 

— sin ( 
r 


■ 9 + - 6 $ 

r 


( 12 ) 


where & includes the non-central terms of F (those that not depend only on the distance 
r and are not directed along the r direction). 



X 


Figure 1. Position vector (r) in the spherical coordinate system used in this paper. 


From ( |TT| ) and (12), we can write the Lagrangian 
& = T-V 

1 -2,1 2d2 i 1 2 • 2 n 12 

= -fir + -fir t) + -fir sin t) <ft — 


1 q 


47ren r ’ 


(13) 


which we call ‘incomplete’ because it contains only the term in (12) that is derivable 
from a scalar (Coulombian) potential. 

As not all the forces acting on the system are derivable from a scalar potential, 
then Lagrange’s equations can be written in the inhomogeneous form [TJ sec. 1.5] 


d 

dt 


d% 

deg 


d& 
dqi 


= Qi 


(t = l, 2, 3, 


(14) 
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where Qi are the generalized forces defined by 

dr 

OQi 


(15) 


Substituting (13), (12) and (15) into (14), we obtain the equations of motion as 


/if — pr9 2 — pr sin 2 9 (ft 2 — 




47T60 r 2 


= 0 


ho 


nr 9 + 2 nrr9 — nr sin 9 cos 9 <f> = - k sin 9 (f> 

47T 


(16) 


ho 


nr sin 9 </> + 2 nr sin 9 cos 9 9(f) + 2 nrr sin 9 <f> = —n sin 9 9. 

47T 


It is worthy of note, however, that the equations of motion (16) could also have 
been obtained from a Lagrangian of the form 


C/> 1-2,1 2q2 i 1 2 • 2 a 12 1 9 , hO 

Jz = -nr H—nr 9 H —nr sm 9 o -1-ncosl 

2 2 2 47re 0 r 47 t 


(17) 


which we call the ‘minimal’ Lagrangian as it is the simpler one that furnishes 


the equations (16) that describe the classical interaction of two dyons without the 


appearance of Dirac strings, according to our generalization of Sokolov’s procedure. 


Now, it is important to note that, since the interaction force (12) is not central, we 


should not expect the vector mechanical angular momentum L to be conserved. As a 
matter of fact, from its expression in spherical coordinates, using (|TTj) again, 

L = pr x v 

= pr 2 sin 9 <f 6 + iir 2 9 <p, 

we can evaluate its temporal derivative by remembering that 


(18) 


— = 96 + sin 9(f 
at 

d 0 • 

— = — 9r + cos( 
df 


0 


d(f> 
d t 


= — sm i 


b <f> 

cos 9(f) <f>, 


(19) 


obtaining traightforwardly 

(£).=»■ 


( —— ) = —nr 2 sin 9 <f) — 2 nr 2 cos 9 \ 

\ dt Je 


2prr sin 9(f> 2 , and 


( 20 ) 


( ) = nr 2 9 + 2 nrr9 — nr 2 sin 9 cos 9 <p 2 

\ dt J6 


Comparing these results with the left-hand sides of the two last equations of 


motion (16) above, we conclude that 


=-^k9 and 


df 


fdL 

\ d t 


47r 


■ n 1 

— k sm 9 d) 
47r 
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or 


d L hq dr 

dt 47t ' dt 


and that the vector 


J = L + j -nr 
47r 


( 21 ) 


( 22 ) 


is conserved. 


The vector J, defined by (22), is known as the Poincare integral of motion, as it was 
first found by Poincare in his previously mentioned work [6|. It can be interpreted as 
the ‘total’ angular momentum of the system because the second term in the right-hand 
side of (22) is the angular momentum of the electromagnetic field, as demonstrated by 
Thomson | G p. 532], 

On the other hand, from the definition (22) and the fact that L is perpendicular 
to r, we also obtain 


J 2 = L 2 + 


ho V 
47T/ 


K 


(23) 


and, since J is conserved and n is a constant, we conclude that the module L of the 
angular momentum is conserved even though the vector L is not. 

We can now use the remaining equation of motion (16) to obtain the conservation 
of the total energy of the system. To do so, it can be rewritten, in terms of L given 
by (18), as 


L 2 

pr - 3 

fir 6 


1 q 


= 0. 


47re 0 r 2 

Now, to proceed further, we may do the trick [I] p. 74] of rewriting its right-hand 
side as a derivative in r and multiplying both sides by r as 

d 
dr 


firr = 


L 2 1 

+ 


9 , 

- r, 


2pr 2 ' 47reor ; 

from what, remembering that d/(r)/dt = (dr/dt)d/(r)/dr = rd/(r)/dr and that 
dr 2 /dt = 2 rr, it follows that 


d 

dt 


h -2 

~ r dt 


+ 


1 q 


2 fir 2 4neo r 


which expresses the conservation of the total energy of the system 


E = —r 2 + 


+ 


q 


2fir 2 47reo r 


(24) 


Notice, now, that the definitions (22) of the J vector and (18) of the angular 
momentum L lead to the result 

ho 


J f — yfir x v + ■ f 

ho 

= K 
47T 

where we used the fact that r x v ■ r = 0. 


( 25 ) 
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If we interpret J r as the projection of J on the direction of r, we can define a as 
the angle formed by r and J given by 


J ■ r 


cos a 


J 


which, as we see from (25) and (22), is constant, with value 
a = arccos(/r 0 K/47r J) 
or, by using the trigonometric identity 


tan(arccos x) = 




x z 


X 


we get from (26) and (23) 


a = arctan 


= arctan 


\/l ~ (/iQ/t/47rJ) 2 
Pok/AkJ 

f 47tJa/1 — (pqk/AttJ) 2 




IMS K 


(26) 


= arctan { 

\ ho« J 

= arctan(47rV / L 2 //Xoft) 

= arctan(47rL//r 0 K). (27) 

Now, being J a vector fixed in the space and the angle it forms with r constant, it 
implies that the motion is limited to the surface of a cone (the Poincare cone ) (Figure[2]) 
of constant half-aperture angle a given by 

a = arccos(/r 0 |k|/ 47rJ), (28) 


with J being coincident with the interior axis if k is positive and with its exterior axis 
if k is negative. A similar conclusion was obtained by Poincare [Bj for the movement of 
an electric charge in the held of a pole of a magnet (equivalent to a Dirac’s monopole) 
and Appel [T6]J for the movement of a electric charge in the held of a magnetic and 
electric pole simultaneously (equivalent to a Schwinger’s dyon). 

The fact that the motion is limited to the surface of a cone of constant half-aperture 
angle a allows us to choose a new spherical coordinate system (r, a, /3), in which the 
vector J coincides with the polar axis, so as to have only the two degrees of freedom 
radial distance (r) and azimuthal angle (/3). 

With this coordinate system, a — 0, while the Lagrangian © reduces to 

= -pr 2 -|—/nr 2 sin 2 a /3 2 -- + ^-ncosa (3 (29) 

2 2 47reo r 47 t 


and the equations of motion (16) result 

1 


pr — pr sin a /3 — 




47reo r 2 


= 0 


ho 


nr 2 sin a cos a (3 = —k sin a (3 
47T 

2 prf sin 2 a /3 2 + pr 2 sin 2 a j .3 = 0. 


( 30 ) 
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In the same way, the angular momentum (18) is now expressed by 
L = —pr 2 sin a /3a, 

the conserved module of the angular momentum as 
L = pr 2 sin a /3, 


(31) 

(32) 

2 ’ 2 pr 2 ' 47reo r ^ 

Now, to arrive at the equation of the orbit, we will follow Goldstein’s procedure jU 
sec. 3.7]. 


and the conserved energy (24) results 

;2 t 2 i 

+ 


pr“ L 2 

£ = — + 


q 


To start with, we can solve (33) for r and get 


r — \ — { E — 


1 q 


L 2 


p 


or, being r = dr/dt, 


, /2 E 

d t = ( — 

p 


47re 0 r 2 pr 2 


Q 


p 2 r 2 


- 1/2 


dr. 


(34) 


(35) 


27reo P r 

For the equation of the orbit, we need the dependence of r upon 6 eliminating the 
parameter t. This elimination can be done by seeing (32) as a relation between d/3 
and d t, in the same way as we did for r: 

Ldt = pr 2 sin a d/3 


or 


d/3 = 


pr z sm a 


-d t. 


The substitution of (35J) into ( |87[ ) yields 
L 


d^ = 


(2E 


q 


pr 2 sin a \ p 27T6o pr 


p 2 r 2 


- 1/2 


dr. 


Now, integrating (38) after slight rearrangements, we obtain 

- 1/2 


P = 


sm a 


2 pE 


pq 


- 4 ) %+t>. 


(36) 


(37) 


(38) 


(39) 


L 2 27T6o L 2 r r^ J r- 

where /3' is a constant of integration determined by the initial conditions and not 
necessarily being the same as the initial angle /3o at time t — 0. 

Finally, changing the variable of integration to u — 1/r, we obtains 


/3 = /3' 


pq 


sm a 


/2 pE 
yT 2 ' ~ 27T6Q L 2 


u 


u 


- 1/2 

2 ' d u. 


This indehnite integral is of the standard form d P- 93] 
dx 1 ( b + 2cu' 


\/a + bu + cu 2 


arccos 




(40) 


(41) 
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A = b 2 — 4 ac, 
where we identify 

a = 2 p,E/L 2 , 
b = —pq/27re 0 L 2 , and 

c = —1. 


Applying (43) to (42), we obtain 


A = 6 2 (1 — 4 ac/b 2 ) 
2 


m 


2ireoL 2 

m 

2nenL 2 


1 + 4- 


2 fiE f 27reo-h 2 

U 2 


2 \ 2 ' 


m 

T 


2 r x 2(47re 0 ) 2 AL 


m 


while applying (43) to b + 2cx results 
b + 2cx = 6(1 + 2 c/bx) 

2n CnL 2 


1 - 2 - 


2t: CnL 2 


Ml 


-u 


Now inserting these results into (|40), we obtains 

'47reo£ 2 w 


P = P - - 


arccos 


sm a 


- 1 


Ml 


1 + 


2(Att€o) 2 EL 2 


Ml 


which we can solve for u = 1/r, obtaining the equation of the orbit as 

Ml 


1 

r 


AnenL 2 


1 + \ 1 + 


32v r 2 dEL 2 


m 


cos [sin a (/3 — f3')\ 
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(42) 


(43) 


(44) 


(45) 


(46) 


(47) 


where we now identify /3 1 as one of the turning angles of the orbit. 

One sees that, except for the sin a term, (47) is very similar to the equation for the 
Kepler problem of the planetary orbits |U p. 93]. As a matter of fact, if the cone has 

(a = 7r/2), the equ; 

[1 + e cos (/3 - fd ')\, 


been degenerated to a plane (a = 7 t/ 2), the equation of the orbit (47) results 

\ = 

r 


47 t er\L 2 


(48) 


which represents a conic curve with eccentricity 

£ = + 32vr 2 elEL 2 //iq 2 . (49) 

Now, remembering that the distance between two points that are differentially 
separated on the surface of a cone with a half angle a is 

(ds) 2 = (dr) 2 + r 2 (sina d/d) 2 , (50) 

while the distance between two points on the plane is given by 

(ds) 2 = (dr) 2 + r 2 (d0) 2 , 


( 51 ) 
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Figure 2. The orbit of the dyon confined to the surface of the Poincare cone with 
half-aperture angle a, according to 


we interpret (47) as representing a conic-shaped orbit confined to the surface of the 
Poincare cone, as shown in Figure |2j 

As in Kepler problem, such conic-shaped orbits divide themselves into bound states 
with elliptic-like orbits (e < 1) and scatterings with hyperbolic-like orbits (e > 1) [T, 
p. 94], We will explore both cases in the next two sections. 
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Let us consider first the bound states (e < 1) of the dyon-dyon system. 


From (49), this corresponds to 

y/l + 32 tt 2 e 2 0 EL 2 /pq 2 < 1 


or 


1 + 32t T 2 e 2 0 EL 2 /pq 2 < 1 
32n 2 e^E L 2 / pq 2 < 0 

which, being everything else positive, implies 
E < 0, 

as expected. 


(52) 


Now, from (33), (52) implies 


nr 2 L 2 1 q 

—~ -b —- 2 + --< 0 

2 2pr~ 47reo r 

1 q pr 2 L 2 

47re 0 r 2 2 pr 2 

which again, being everything else positive, implies 

q< 0 , 

i.e., an attractive interaction with 

1 Igj »r 2 L 2 

47re 0 r 2 2 pr 2 

Being elliptic-like orbits, we can calculate, from (47), its two turning points (points 
of least or greatest distance of one dyon from the other), r min and r max , as 


(53) 


T min. — 


Vmnrr = 


47reoL 2 

v\q\ 

4ne 0 L 2 

v\q\ 


(1 + e) 


(1-e) 


(54) 


from which we can obtain the semi-major axis a and semi-minor axis b as 

Vmin “b r rnax 
a =- 


b = \fa 2 { 1 — e 2 ). 

However, it is more useful to have a and b expressions in terms of the energy of the 
system. To start with, we make use of the fact that the radial velocity r is zero at those 


turning points. Therefore, from (33), the conserved energy at those points becomes 

L 2 


E = 


1 q 


2 pr 2 4ne 0 r 
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which can be rewritten in the form of a quadratic equation as 
1 a L 2 


r 2 — 


r — 


= 0 


47reo E 2pE 

having the turning points r min and r max as its roots. Now, it is well known that the sum 
of the roots of a quadratic equation equals the negative of the coefficient of its linear 
term. Therefore, 


T T'max 


q 


47T60 E ’ 


and we obtain the semi-major axis a as 
Tmin T r max 1 Q 


a = 


Sneo E 


(55) 


and the semi-minor axis b as 
b = a/ci 2 (1 — e 2 ) 



1-1 + 


32? r 2 e 2 0 EL 2 ^ 


ME 


\QA'K 2 e^ ) E 2 / 
L 


/ 32t r 2 e 2 0 EL 2 \ 
V M 2 ) 


(56) 


y/tyAE | 

Notice that the existence of r mm and r max does not necessarily mean that the orbit 
is closed but only that it is ’bounded’ by those limiting distances. 

To analyse the closedness of the orbits, we have to see that the term 


cos [sin a {(5 — (5’)] in (47) implies that the ratio between the periodicities of the radial 
coordinate (r) and the azimuthal angle (/3) is given by sin a. Consequently, it is this 
parameter that will determine if the orbit is closed or not. 

To understand that, we have to consider that, at each revolution, the dyon describes 
a portion sin a of the ellipse. Now, the orbit will be closed if, say, after n revolutions, it 
will have completed exactly m ellipses, that is, if and only if sin a is rational 

(57) 


m 

sin a = —, 
n 


with m and n relatively prime numbers and m < n HI p. 91]. 

Furthermore, the numbers m and n define the topology of the orbits in terms of 
number of double points [H, p. 28], classifying them into families of the same topology. 


5. Dyon-dyon scattering 

After the dyon-dyon bound states, seen in the previous section, let us now study the 
classical scattering (e > 1) of one dyon by another. 

For this study, we need to obtain the so-called cross-section for scattering in a given 
direction cr(fl) and to arrive at it we will follow Goldstein’s procedure [U sec. 3.10]. 
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to 


From (49) and a reasoning similar to the one that leaded to (52), this corresponds 
E> 0, (58) 


which, from (33), analogously to (53), implies a repulsive interaction, i.e. 

q > 0. (59) 

To start with, let us define Uo as the velocity of the dyon when it is at an infinite 
distance (r —y oo) from the other. Then, from (33), the conserved energy reduces to the 
kinetic energy at that point: 

1 


e = yy. 


( 60 ) 


We can now define d as the distance of closest approach (periapsis) and impose the 
following condition on fd' 

/3' = p( r = d) = 0 (61) 


into the general equation of the orbit (47), obtaining 
1 
r 


pq 


47t enL 2 


[1 + ecos (fd sin a)] 


( 62 ) 


( 63 ) 


with e now given by substituting (60) into (49) as 
£ = \Jl + 32n 2 el(in%/2)L 2 /pq 2 

= yj 1 + ( 4 tt e 0 ) 2 L 2 v$/q 2 , 

a result similar to that obtained by Schwinger m- 

We can now study the dyon-dyon scattering by defining the scattering angle 0, the 
angle between incidence f t and scattering f/ directions at an infinite distance from the 
origin, as 

cos© = —Vi ■ Vf. (64) 

Expressing r in the coordinate system (r, a, (3) by 

r = sin a cos fd i + sin a sin fd j + cos a k (65) 


with a constant and given by (28), we obtain, from (64), 

cos© = —[sin 2 a (cos fdi cos/3/ + sindj sin/3/) + cos 2 a], 

where fdi and fdf are the azimuthal angles of the directions of incidence and scattering 
at an infinite distance from the origin, respectively, and, by making use of a few 
trigonometric identities such as 

cos 2 6 = 1 — sin 2 9 (66a) 

cos($! — 6*2 ) = cos^x cos# 2 + sin 6\ cos #2 ( 666 ) 


1 — cos 9 


sin 2 (#/2) = 

z 

cos(#/2) = ±y/ (1 + cos#)/2 


(66c) 
(66 d) 
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cos© = —[sin 2 a cos(/% — /3f) + (1 — sin 2 a)] by (66a) and (66© 

= sin 2 a [1 — cos(/3; — /?/)] — 1 
= sin 2 a {2sin 2 [(/3j — /3/)/2] j — 1 by (66c) 

= 2 {sina sin[(ft - /3 / )/2]} 2 - 1, 
and, finally, by (66d), 

cos(0/2) = sin a [sin [(/3j — (3f)/ 2] | . (67) 

To proceed further, we need to calculate /© and /3/, the azimuthal angles of the 
directions of incidence and scattering at an infinite distance from the origin. Notice 
that (1 /r) r . 


= 0 and, therefore, from (|62|), we get 
fiq 


[1 + e cos (/3 sin a)] 


= 0 


47re 0 T 2 

and we can evaluate cos(/3 sin a) r _s.oo as 
cos(/S sin a),—s-oo = —1/£ 

= -[l + (4 ^„) 2 LW]- 1/2 

= — [l + (47rei)) 2 (/io^ tail a/4ir) 2 i;o/5 2 ] 
= — {l + (k 2 Hq/c 4 ^ 2 ) tan 2 a) 1,/2 


( 68 ) 


where we used (63), (27), and = 1/c 2 . 

Now, by making use of the trigonometric identity 

cos(arctan:r) = 1 


vTT 


x 


2 ’ 


into (68), we obtain 

(j3 sin ck), 


= arctan [{kvq/ c 2 \q\) tana] 
= arctan (fctana), 

where we introduced the parameter k defined as 
k = KV 0 /c 2 \q\. 


Therefore, from (67), we finally obtain the scattering angle 0 as 
cos(0/2) = sin a |sin($/ sina)| , 

where, for compactness, we introduced the parameter <f> defined by 
<f> = (/3 sin a^r-Kjo = arctan (k tan a), 


(69) 

(70) 

(71) 


from the result (69). 

We present in Figure [3]graphics of the functional dependence of 0 upon a for various 
values of the parameter k. In this figure, it is also displayed a graph for monopole- 
electron scattering, obtained by making e 2 = g\ — 0 (which, from ([9]), results in q — 0 
and k = —eg/c ) in (69). 

















Figure 3. Graphs of the functional dependence of the scattering angle 0 upon the 
half-aperture angle of the cone, a , according to for the monopole-electron pair 
and for various values of the parameter k = nvo/c 2 q. 


As the scattering orbit is asymptotic to r* and ff, which define a plane, we can 
calculate the elastic differential cross-section [U p. 108] for the dyon-dyon scattering, in 
terms of the impact parameter s (the perpendicular distance between the centre of force 
and the incident velocity) [IJ p. 107], by the equation 

(72) 


da 

sds 

- V s 

d(cos 0) 

dO 

d(cos 0) 

s 

ds 


since, as seen from the graphs in Figure [3j in general, a is a multivalued function of 0 [Q~ 

p. 111]. 

As Vo is the incident velocity and perpendicular to the impact parameter s, we can 
express the angular momentum L, from (18), as L = psv o and also substitute (27), 
resulting 


s — L/ (fjivo) 

= (/iok/47t/zu 0 ) tana, (73) 

Now, using d/(a)/ds = (da/ds)d/(a)/da = (ds/da) _1 d/(a)/da and tan' 6 = 
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seer 9 = cos 2 9, we can rewrite (73) as a relation between ds and da as 

-i i 


f JM o«_\ 


ds [\i7iiw 0 J cos 2 aJ da’ 


(74) 


and use this and (|73j) to evaluate the derivatives with respect to s in (|72j), obtaining the 

1 d(cos@) 


elastic differential cross-section as 


da 

dn 


E 


f \ 

\47r/ino J cos 2 a 


E ( 9-qK \ 
l Airline I 


sma 


\47rpu 0 J cos 3 a 


da 
d(cos 0) 


-l 


( V 

- tana 

\47r/in 0 / 


da 


-l 


( k0K N 


sin(2a) da 

V 4? rpu 0 , 

/ 2 cos 4 a 

sin 0 d@ 


(75) 


where, in the last step, we used the trigonometric identity sin(2a) = 2 sin a cos a. 

Particularly interesting are the cases where the parameter k dehned in (70) is a 
multiple of tt. One can see that by expanding (72) with (69) around a = 0, using the 
approximations tan 0 = sin 6 and arc tan 0 = 6, as 

$ = arctan [{kvq/ c 2 \q\) tana] 

= (KVo/c 2 \q\) sina, (76) 


and using this result in (71), obtaining 

cos(0/2) = sina |sin{[(«;no/c 2 |g|) sina]/sina}] 

= a |sin(«n 0 /c 2 |g|)| , 

from which, by means of the trigonometric identity cos(20) = 2cos 2 0 — 1, we obtain 

cos(O) = 2a 2 sin 2 (/cu 0 /c 2 |g|) — 1. (77) 

With these approximations and the usual sin 6 = 6 and cos 6 = 1 ones, the cross- 
section (75) results 

da a _>o 


a —>0 

—y 


a^O 

—y 


E f p 0 K \ 

\ 47rp.no/ 01 

Ot N ' 

( Hq K V 


d[2a 2 sin 2 (s:no/c 2 |g|) — 1] 


da 


a sin 2 (KVo/c 2 |g|) 


-l 


i-i 


(78) 


\8nnvoJ sin 2 (/vu 0 /c 2 |g|) 
and, therefore, when a —> 0, a second order pole occur in the cross-section whenever 

k — rnr (rz = 1,2,3,---). (79) 


To evaluate the differential cross-section (75), we found it convenient to introduce 
the new variable 


( = 2<f>/ sina, 


(80) 
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where $ is given by (72) and (69) 


Now, from (71) and again by means of the trigonometric identity cos(20) = 


2 cos 2 6 — 1, we have 

cos(0) = 2 sin a 2 sin 2 (£/2) — 1. 


(81) 


From this result, we can evaluate the derivative on the right-hand side of (75) as 

-l 


1 dct 
sin 0 d@ 


d cos 0 
da 

= sin(2a) 


(1 — cos () + sin (tan a 


d(C/2) 

da 


Then, from ( [80] ) and (72), we can evaluate the derivative of (/2 as 
d(C/2) kv o cos 2 [(C/2) sin a] 


da 


q sm a 


cos 2 a 


— (C/2)cot a 


and obtain, from (75), the differential cross-section as 
/ no k \ 2 


da 


where 


9(0 = 


\47T/it’o J 
1 


9(C), 


cos 4 a |2(1 — cosC) — C s i n C + (sin. C/ sin a) sin(C sin a)| ’ 


(82) 

(83) 

(84) 

(85) 


with C given by 

In Figure [4] we present graphs of the differential cross-section for various values 
of the parameter k. In this figure, we also show a graph for the electron-monopole 


scattering obtained, again, by making q = 0 and k = —eg in (85) 


One observes from Figure [4] that the cross-section becomes infinite for some values 


of 0 and falls abruptly thereafter. From (75), we see that it happens, besides 0 = 0, in 
one of the following cases 


d0/da = 0 

0 = 7T but sin(2a) ^ 0. 


( 86 ) 


This phenomenon is very similar to what occurs in the optical scattering of sunlight 
by raindrops and, because of this similarity, these two conditions are referred to as 
rainbow scattering [H p. Ill] and glory scattering [TJ p. 114], respectively. 

In the case of rainbow scattering, we have, from 


tan 






+ 


sin(2$ r ) 


= 0, 


(87) 


k sma r y sma r sin(2a r ) cos a r 
where the suffix r was appended to indicate that it refers to the rainbow scattering. 
Equation (87) can be solved numerically. In Table [lj we present rainbow a r angles for 
two values of the k parameter defined in ( jT0| ). 

For the glory scattering, we have, from ( |7Tj ), 

<3> g /sin a g = wr ( 77 . = 1,2,3,---), 


= nn 


( 88 ) 
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1.5 

1.0 

0.5 

0.0 




Figure 4. Differential cross-section g(Q) scattering, according to ( |85| , for the 
monopole-electron pair and for various values of the parameter k = nvo/c 2 \q\. 


Table 1. Rainbow angles a r (d0/da = 0) obtained from (871. 



k 

a T 

2tt 

0.230 

3tt 

0.104 


0.280 


and, from (|72j), the condition 

tan(mrsina g ) = ktana g (n = 1,2,3,-••), (89) 

which can be solved numerically for the desired value of the parameter k. In Table [2j 
we present glory a g angles for two values of the parameter k. 

Let us now consider what happens at small-angle scattering (a —> vr/2). From (71) 
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Table 2. Glory angles a g (0 = n) 


k 

a g 

2tt 

0.394 

37r 

0.156 


0.445 


obtained from (89). 


and using the trigonometric identity sin arctan 6 = x/y/l + x 2 = 1 /\/l + x ~ 2 , we have 


cos(@/2) = sin a sin [arctan (fc tan a) 
sin a 


a/1 + (k tan a) -2 


(90) 


Now, using (27) and expressing the parameter k in terms of the impact parameter s 
by means of (73), we have 

tana = (47 t//xo k )£ 

= ( 47 T /hqk){hvqs) 

= (4:TTfXV 0 /fl 0 n)s, (91) 

and, by means of the trigonometric identity sin 0 = 1/y/l + cot 2 6, we obtain 

sin a = l/\/l + (/io^/4vr/in 0 s) 2 . (92) 


On the other hand, from (70) and (91), we get 
A; tana = {kvq/ c 2 \q\){A'K pvo/ iiqk)s 
= (47re 0 /iWo s /M), 

where we used po^o = 1/c 2 once more, and, therefore, 

1 1 


y/l + (A: tan a) -2 y/l + (g/47re 0 /ifo s ) 2 


(93) 


(94) 


Substituting (92) and (94) into (90), using the approximation 1/y/l + x = 1 — x/2, 

1 


and keeping terms only up to first order, we obtain 

1 


cos(©/2) = 


\J 1 + (g/47re 0 /iWo s ) 2 \A + (noK/^nfiVosy 


1 

1 - - 




2 \47reo^UoS 


= 1 - - 
2 


Q 

47re 0 /ir'o s 


+ 


= 1 - - 


2 (47re 0 /ino) 2s2 


1 _ 1 / /JqK 
2 \47Tyuno-s 


f 

\4niiv 0 s 


K 


— 1 + ( - 


(95) 
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Now, considering the approximation cos(0/2) = 1 — (@/2) 2 /2 and comparing this 


expression with (95), we have 


( 0 / 2) 2 = 


1 


(47re 0 /ino) 2 s 2 


or 


s 2 2* 


(0/2) 2 (47re 0 /mo) 2 


+ 1A 2 

v 0 


hi 

Vo / + Vc 2 


(96) 


(97) 


Thus, from (72) and the approximation sin 9 = 9 , we can compute the differential 


cross-section of the dyon-dyon scattering at small angles as 


dcr 

dn 


= s 


ds 


d(cos 0) 


= s 


= s 


ds 


d(s 2 ) 

d[(©/2) 2 ] 

d(0/2) 

d(s 2 ) 

d[(©/2) 2 ] 


d(©/2) 

d(©) 

d( S 2 ) 

-l 

d[(e/ 2 ) : 

] 

d(©/2) 

d(cos0) 

ds 


d(0/2) 


d(0) 

d(©) 


d(0) 


d(cos0) 

1 d(s 2 


d[(©/2) 2 


= s|2s| 1 |2(0/2)| - | —sin0| 


-l 


d(s 2 


d[(©/2) 2 


4 | sin @| (47reo/iUo) 5 


- +(4 

v 0 ) Vc 2 


d[(0/2) 2 ] \ (0/2) 


4 (47r6 0 pv 0 y 


4 (47re 0 /iUo) 2 


— + Hr 


— + Hr 


(0/2) 4 

csc 4 (0/2) 


(98) 


which is a generalization of the Rutherford formula for the scattering of a particles by 
atomic nuclei HI p. 110]. 

A relativistic extension of the entire calculation previously done is possible, in a 
simple way m cap. 4]. The importance of a relativistic extension for the treatment 
of dyon-dyon system stems from the large value of the coupling constant for magnetic 
charges (see (J3])) . In the same way as in the non-relativistic case, the relativistic dyon 
is confined to the surface of a cone of half-aperture angle given by (27). The bound 


states correspond to conic-shaped orbits confined to the surface of the Poincare cone 
(Figure [2]), with a similar condition for the orbit closing. For the relativistic scattering 
case, the orbits are hyperbolas confined to the surface of the Poincare cone and this 
system exhibits glory and rainbow scattering. 

It is even possible to apply the well-known Sommerfeld (semi-classical) quantization 
rule Ha p. 283] to the dyon-dyon system. By doing it, we have shown the analogy 
between this system and the hydrogen atom and that the former can be considered 
its generalization [2] cap. 3]. Furthermore, we obtained, for both the non-relativistic 
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and relativistic cases, an energy spectrum that reasonably approximates the quantum 
spectrum obtained by Pereira [19], 


6. Concluding remarks 

In this work, thanks to a Lagrangian without Dirac strings , it was possible to make 
a classical study of the dyon-dyon system, with the electron-monopole system as a 
particular case. 

The orbit equations were interpreted as representing conic-shaped orbits confined 
to the surface of the Poincare cone. While the electron does not form bound states with 
the magnetic monopole, the dyon-dyon system exhibits non-planar stable elliptic-like 
orbits, and we presented the conditions for them to be closed. 

Through the orbit equations, it was also possible to study the classical scattering for 
these systems. We showed that the elastic differential cross-section shows divergences 
that are usually denominated glory and rainbow scatterings by their similarities with 
Optics. We also showed that the differential cross-section of the dyon-dyon scattering 
at small angles is a generalization of the Rutherford formula. 

Finally, we want to stress that a relativistic extension of this study can be done in 
a simple way, as well as a semi-classical quantization via Sommerfcld rule. 

We hope that these results may arouse students’ and teachers’ interest and 
contribute to Classical Mechanics courses. 
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